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Recent discussions of RHIC data emphasized the exciting possibility that the matter pro- 
duced in nucleus-nucleus collisions shows properties of a near-perfect fluid. Here, we aim 
at delineating the applicability of fluid dynamics, which is needed to quantify the size of 
corresponding dissipative effects. We start from the equations for dissipative fluid dynam- 
■ ics, which we derive from kinetic theory up to second order (Israel- Stewart theory) in a 

\ systematic gradient expansion. In model studies, we then establish that for too early initial- 

ization of the hydrodynamic evolution (tq<1 fm/c) or for too high transverse momentum 
(pt^I GeV) in the final state, the expected dissipative corrections are too large for a fluid 
description to be reliable. Moreover, viscosity-induced modifications of hadronic transverse 
momentum spectra can be accommodated to a significant degree in an ideal fluid descrip- 
tion by modifications of the decoupling stage. We argue that these conclusions, drawn from 
model studies, can also be expected to arise in significantly more complex, realistic fluid 
^vq ■ dynamics simulations of heavy ion collisions. 
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Why is it interesting to characterize dissipative effects of the dense QCD matter produced in 
ultra-relativistic nucleus-nucleus collisions at RHIC or at the LHC? The experimental heavy ion 
programs aim at establishing properties of QCD matter at the highest energy densities attained in 
the laboratory 0, 0, 0, Q|. Shear viscosity, bulk viscosity, heat conductivity or the conductivities of 
conserved charges are such properties of QCD matter. If unambiguously extracted from data, they 
are prime candidates for the next compilation of the Particle Data Group. They are of fundamental 
interest, since they are not mere material constants, but they are - at least in principle - computable 
from first principles @, @, 0, S 0| • 

How can one extract dissipative transport coefficients from data? For this, it is a prerequisite to 
have a dynamical theory of nucleus-nucleus collisions, which can be compared to data and depends 
on dissipative properties. Dissipative hydrodynamics is this theory. It describes liquids, which 
deviate locally from a fully equilibrated, 'perfect' one. Dissipative hydrodynamics is not applicable 
to general non-equilibrium evolution. Deviations from ideal fluid dynamics must be sufficiently 
small for the gradient expansion underlying dissipative hydrodynamics to be valid. On the other 
hand, dissipative effects must be sufficiently large to be measurable. These two requirements limit 
the applicability of dissipative hydrodynamics to heavy ion collisions. They will complicate any 
attempt to determine transport coefficients from data. 

Remarkably, simulations of Au-Au collisions at RHIC based on ideal fluid dynamics [lfl . Hll . 
12l Hi . 3] reproduce the observed large size and impact parameter dependence of elliptic flow for 
sufficiently central collisions (b < 7 fin) near mid-rapidity 15, 16, 17]. They also account for the 
transverse momentum dependence of hadronic spectra up to pj- < 1.5 GeV, and they reproduce 
the gross features of the particle species dependence of these spectra 2 3- However, these ideal 
hydrodynamic descriptions of RHIC data require a very short thermalization time, tq < 1 fm/c |2C i. 
for which multiple scattering models cannot account naturally in the weak coupling regime 21, 22]. 
While this adds to the paradigm that high-energy heavy ion collisions create a strongly coupled 
plasma ijj 3 13- 2i3- 2?], 13] , one cannot exclude a priori the existence of collective mechanisms 
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(such as plasma instabilities 0, EJEil HIS), 

which may account for fast equilibration in the 

weakly coupled regime. 

Ideal hydrodynamic simulations of Au-Au collisions at RHIC also gave support to the conclu- 
sion jj], H, El, [ij that for realistic initial conditions, the observed elliptic flow exhausts the theoretical 
upper limit predicted by ideal fluid dynamics, thus indicating that a perfect liquid with negligible 
dissipation has been created at RHIC |23L Efl|. However, this statement must be qualified, since 
recent work j3] has identified a class of a priori realistic initial conditions, for which ideal fluid dy- 
namics significantly overestimates the elliptic flow measured at RHIC. Also, it has been questioned 

I L 

that the observed elliptic flow is indicative of full equilibration 35]. This illustrates the major 
limitation of hydrodynamic concepts to heavy ion collisions: while the dynamical description is 
parameter free, rather precise knowledge of the initial conditions is essential for the predictive or 
interpretational power of the theory. This is so for ideal fluid dynamics already. But since the ap- 
plicability of dissipative fluid dynamics depends on precise knowledge about the closeness to local 
equilibrium, one may expect that control over initial conditions becomes an even more sensitive 
issue for the characterization of dissipative properties. 

The main motivation of this study is to delineate the region of applicability of (ideal as well 
as dissipative) fluid dynamics in heavy ion collisions. To this end, we consider dissipative hydro- 
dynamics up to second order in the gradient expansion around ideal hydrodynamics, the so-called 
Israel-Stewart theory [3(|. This theory has been applied to viscous hydrodynamics only recently 
0, EI (see however 0, subjected to further approximations which we discuss. The 
present work follows the standard perturbative logic, that to check the validity of the zeroth or- 
der (i.e. ideal hydrodynamics), one should check the size of higher orders in the perturbative 
expansion (here: gradient expansion). An expansion to second order is necessary, since first order 
dissipative hydrodynamics implicitly assumes vanishing relaxation times and is known to show 
therefore acausal artifacts |4^|. Beyond quantifying the 'closeness' of the time evolution to ideal 
hydrodynamics, the motivation for studying the gradient expansion of dissipative fluid dynamics is 
of course to gain access to the fundamental dissipative properties of the produced matter, namely 
its transport coefficients. 

Our work is organized as follows: In Section [Hi we specify the equations of motion of second or- 
der dissipative hydrodynamics. We have found these equations in the literature subjected to further 
sometimes ad hoc approximations. In Appendix we rederive these equations from Boltzmann 
transport theory. This allows us to discuss explicitly differences between existing formulations. 
We then specify to a simple model inspired by heavy ion collision in Section II I H and we discuss 
analytical results for this model. In Section ITvl we numerically solve the second order dissipative 
hydrodynamic equations. This leads to a lower bound on the hydrodynamic initialization time tq 
and upper bound pj. lt for the transverse momentum up to which a fluid description is reliable. We 
then proceed to numerically evaluate hadronic transverse momentum spectra via the Cooper-Frye 
formalism (4^, both for ideal and second-order dissipative hydrodynamics. These results illustrate 
the extent to which particle spectra emitted from a dissipative fluid may be mimicked by an ideal 
fluid dynamics description. We discuss our main conclusions in Section 



II. DISSIPATIVE HYDRODYNAMICS 

From the Boltzmann transport equations, one can derive the hydrodynamic equations of motion 
including dissipative corrections to second order (see Appendix 

(e + p)Du^ = V^p - A^VJT* 7 + W v Bu v , (1) 
De = -(e + p)V^ + ^n^(V^), (2) 
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m A^DU a(3 + W u = r?( W) - 2T U U a ^co v l . (3) 

Here, the energy density e and the pressure p are related by the equation of state. The vector 
denotes the local fluid velocity, and is the shear tensor. In the limit of vanishing relaxation 
time Tn, the shear tensor is given by instantaneous information about the gradients of the fluid 
velocity, IP 41 ' = rji^J^u 11 ). This is the definition of shear viscosity in first order dissipative fluid 
dynamics. In principle, dissipative corrections depend also on heat flow and bulk viscosity, as well 
as on the corresponding relaxation times. However, the effects of the latter are expected to be 
much smaller on general grounds. To arrive at a transparent discussion, we shall neglect them in 
what follows. We use the following notation 



d^u v 


= d^ + T^u a 


(4) 


D 


= 


(5) 


v M 


= A^ u d u 


(6) 


A** 




(7) 






(8) 




= A^B V + AyB^ — — A. f j il/ A a B a 


(9) 




= i (A^Bv + A u B fl ) 


(10) 


[A^B U ] 


= X - (A^B V - A u B fl ) , 


(11) 



where is the covariant derivative, T^p are the Christoffel symbols and to' 11 ' is the vorticity tensor. 

A. Approximation to hydro dynamic equations 

In recent studies of second order dissipative fluid dynamics one uses for the equation of 
motion of the shear viscosity 

m DIL» u + = rfiW) . (12) 

To understand how this is related to expression (|3J), we expand Eq.(j3J) in the form 

r n [DW + (um u a + u v W a ) (Du a )} + W u - = -2r n n a ^c;^ • (13) 

For 1)12|) and ()13j) to agree, we require that H a ^uj U a = and H^ a Du a = 0. By construction, the 
shear tensor is orthogonal to the fluid velocity, Ii^ a u a = 0, and these two conditions can be written 

as 

u v » = 0, (14) 
Du u = 0. (15) 

Thus, to arrive at Q12|). one has to assume that the liquid is vorticity free. Moreover, one assumes 
ad hoc that Du u vanishes up to second order in a gradient expansion. This second requirement 
is not satisfied for a general flow field, but it holds for the simple example studied in Section IIIII 
below. 
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We note that for a consistent treatment, the assumptions (|14|) and (|15[) have to be applied to 
all three equations of motion (^p)-©- The equations of motion ®-(J3J) are then replaced by Eq. 
(fT2|l and 



D e = -( £ + p)V/ + ~ir^ < > 



(16) 
(17) 



The equations of motion Q-© were obtained from transport theory. An alternative derivation 
starts from the most general second order gradient expansion of the entropy density s. Since entr opy 
cannot decrease, this leads to the requirement that the four-vector of the entropy density [4J| 
satisfies 



TV 



n 



i 

2?7 



(18) 



where = From this relation it is obvious that possible terms in the square bracket orthogonal 
to cannot be found by this "entropy- wise" derivation of the hydrodynamic equations. As we 
have seen above, it is by these terms that Eq. (fT3)) differs from Ea.p2jl. 



III. IDEAL VERSUS VISCOUS FLUID DYNAMICS FOR A SIMPLIFIED MODEL 

For general initial conditions, the hydrodynamic evolution according to Eqs. ©-© is com- 
plicated and requires a (3+l)-dimensional numerical simulation. Here, we consider a simplified 
model, which incorporates several characteristic features of a heavy ion collisions, while lacking 
many of the technical complications of the most general solution. In this Section, we specify the 
geometry, initial conditions, equation of motion, equation of state and viscosity for this simple 
model. 



A. Geometry and initial conditions 

For the discussion of relativistic heavy ion collisions, it is useful to transform to radial coordi- 
nates in the transverse plane, and to proper time r and space-time rapidity ij, 

t = \A 2 — z 2 , r> = arctanh-, r = J x 2 + y 2 , 6 = arctan— . (19) 

t v x 

The metric in these coordinates is diagonal, 

g»v = (9TT,9rr,9<l><f>,9rm) = -1 > _r-2 > ~ T J • ( 20 ) 

The only non- vanishing Christoffel symbols are = r, T^ v = 1/r, = — r, vf^ = 1/r. 

To arrive at a simple model, we assume that the initial conditions are longitudinally boost- 
invariant, i.e., the energy density has a vanishing gradient in r\ and the initial fluid velocity has the 
profile 

u v = { u T ,u T ,u' t, ,uP) = (1,0,0,0). (21) 

The hydrodynamic evolution preserves this boost-invariance, which in ultrarelativistic heavy ion 
collisions is expected to be realized near mid-rapidity j4f|. 
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Second, we assume that the initial conditions show neither flow nor density gradients in the 
transverse direction. This implies that the r- and ^-dependence of the hydrodynamic evolution 
becomes trivial. This condition is realized for matter of r-independent density and infinite trans- 
verse extension. For realistic collision geometry, the edge of the transverse overlap region will 
always show sizeable transverse density gradients, which drive sizeable transverse flow gradients in 
subsequent steps of the dynamical evolution. However, for matter close to the transverse center 
of the collision, approximate r-independence of density and flow can be realized for early times. 
This is also seen by multi-dimensional simulations of ideal fluid dynamics, which reveal for an 
extended time period approximately one-dimensional expansion dynamics. In this sense, we view 
the one-dimensional model resulting from the above-mentioned approximations as a model which 
retains characteristic features of a fully multi-dimensional simulation of heavy ion collisions. 



B. Equation of state and viscosity 

To close the hydrodynamic equations of motion, one has to specify an equation of state. Here, 
we consider the ideal one, 

e = 3p. (22) 

For numerical calculations, we shall use the perturbative result for the pressure of the Yang-Mills 
theory to leading order in the coupling g 2 

tt 2 T 4 2 7N c N f 



p =l—(2(N c -iy + ^+) . (23) 



We consider QCD (N c = 3) with Nf = for simplicity, where 

a = P /T*= 8 -^. (24) 

For numerical calculations, we also need to specify the dimensionless ratio r]/s. The perturbative 
result, given more explicitly in Appendix [UJ is to leading logarithmic accuracy of the form 7] = 

F ln/xVmg - F ° r QCD = ^ and Nf = °' [t takeS the f ° rm ^ 

_ 1 27.126 

r?/S "4a 5 4 ln(2.765 5 - 1 )- ( 5) 

2 

According to this expression, the range a s = §- = (0.2,0.4)) corresponds to rj/s ~ (1.1,0.73). 
However, often the logarithmic correction in Ea. (|25|) is assumed to be ln/j,*/m,£) ~ 0(1)- Then, 
one finds for the same values of the coupling constant the range rj/s ~ (0.61,0.15). This is one of 
many illustrations of the significant numerical uncertainties related to the use of these results. On 
the other hand, there is a conjecture 0, that in all thermal gauge field theories rj/s > ~ 0.08. 
This lower bound is realized in the strong coupling limit of highly super symmetric thermal Yang- 
Mills theories, and it may be of relevance to the strong coupling regime of finite temperature QCD 
near T c . 

We emphasize that we do not advocate the applicability of perturbative results to the dense 
QCD matter produced in nucleus-nucleus collisions. The sole purpose of the above discussion was 
to survey the range and uncertainties of the numerical input, suggested by different perturbative 
and non-perturbative calculational approaches. 
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C. Equations of motion 

In our model, the rj- and r-independence of initial density and flow profile implies strong sim- 
plifications. Making use of (V r ii T ) = and U TT = (following from u a (S7 a u^) = = u a Il a ^), and 
the trace condition 11^ = 0, one finds after some algebra 

in^(V^) = -rn^. (26) 

Since (|14|) and ()15|) are satisfied, the equations of motion take the form 

= Vp-A^n w , (27) 
De = -{e + p)V u u^ -rn w , (28) 
t u DU^ + W = n^vP) . (29) 

Ea. (|28|) can be simplified further by realizing that V pU 11 = 1/r. Moreover, explicit calculation of 
(|2Tj) results in 

V> = . (30) 

Thus, the present model describes matter with vanishing pressure gradients. However, the model 
allows for a non-trivial time-dependent evolution of density since the matter is subject to Bjorken 
boost invariant expansion. 

We now discuss more explicitly, how leading and higher order dissipative effects enter the 
equations of motion. In our model, all information about dissipative effects resides in the component 
n^* 7 of the shear tensor. 



1. Zeroth order - Euler equations: the perfect liquid 

In the absence of dissipative effects, il j?7) = 0, one recovers the equations of motion of an ideal 
fluid 

d T e = - { -^. (31) 
r 

For the equation of state e = 3p, this has the solution 4 



T(t) = T, 



_ . 1/3 
T J 



(32) 



In what follows, this result defines the baseline on top of which dissipative effects have to be 
established. 



2. First order - Navier- Stokes: reheating artifacts at early times 

Dissipative corrections to first order in the gradient expansion are recovered by setting the 
relaxation time rn = 0. This leads to = rj(SJ^u v ), where r\ is the (temperature dependent) 
viscosity. The resulting equation of motion reads 

a e + p 4rj 

d T e = — + ^ • (33) 
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For the equation of state and viscosity defined in Section IIIIBI the solution of this differential 
equation is known analytically for the case of constant rj/s [371 l46j| 



T(t) = To 

' T 



TO X 1/3 



277 

1+ ' 



3st T 

This function has a maximum temperature at time 




(34) 



(1 st T q \- 3/2 
W = 7 tt ^- + -— J . (35) 

For times r > r max , the temperature decreases with time, as expected for matter undergoing 
expansion. For early times r < r max , however, the solution shows an unphysical reheating. We 
note that first order dissipative fluid dynamics is known to show unphysical effects |42| . 

From a pragmatic point of view, one may ask how the numerical importance of this reheating ar- 
tifact depends on the time tq, at which one initializes the dissipative evolution. For one-dimensional 
viscous first order hydrodynamics, the entropy evolves according to d ^T s ' = ^0-. For a perfect 
liquid, entropy does not change. As a consequence, for dissipative corrections to be small, one 
requires that the increase of entropy is small compared to the total entropy in the system, i.e., 
2 -C tT must hold throughout the dynamical evolution. For initial time tq and initial temperature 
To, we thus require 2 -c tq Tq, and we find 

3/2 



T max J- 



—) , for ^«t T . (36) 



VtqTq \s J ' s 



So, if 2 -c ro To, then T max <C to and one can expect that the unphysical reheating effect is not 
seen during the time studied by the evolution. One may then hope that the calculation returns 
reasonable results, although the absence of an obviously unphysical behavior by no means implies 
that. A check by including higher order corrections appears to be advisable even in this range. On 
the other hand, for fixed value ^ < 1, one can always find sufficiently small times To so that the 
entropy density increase per unit time is large initially. For these early times, dissipative effects 
will be large, and will lie outside the validity of the relativistic Navier-Stokes equations. This is 
seen in particular in the appearance of an unphysical reheating effect at early times. 



3. Second Order - Israel- Stewart theory: causal dissipative hydrodynamics 



For the boost-invariant model defined in Section [ill Al the equations of motion to second order 
in dissipative gradients read 37J, |47 ] 



e + P 1* 
- + -$, 

T T 

$ 2 
+ 



tit 3t/3 2 



R - Tn 



Here, we have introduced 



- 2 rf m . 



(37) 
(38) 

(39) 



The ratio of relaxation time over viscosity can be calculated for a massless Boltzmann gas by 
using Ea. (|A24|) and determining the pressure p. One finds 



(40) 
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FIG. 1: Temperature evolution in ideal hydrodynamics (dashed line), first order dissipative hydrodynamics 
(relativistic Navier-Stokes, dash-dotted line) and second order dissipative hydrodynamics (Israel-Stewart, 
full line), rj/s — 0.3. For sufficiently early initialization time To, the relativistic Navier-Stokes equation gives 
rise to an increase in temperature and thus total system energy, which is not seen in the Israel-Stewart 
theory (see text for details). 

For a Bose-Einstein gas, this relation is only approximately true. It is modified by the factor 
~ 1.024 (see Appendix [BJ). The equations have to be solved numerically which we do for the 
equation of state e = 3p. We will scan a range of initial conditions for the temperature and initial 
time To, ro as well as the strength of dissipative effects controlled by rj/s; for simplicity we always 
set $(ro) = in the following. 

IV. NUMERICAL RESULTS FOR SECOND ORDER VISCOUS HYDRODYNAMICS 

We now turn to the numerical solutions of the simplified hydrodynamic model, defined above. 
We shall compare results for ideal fluid dynamics to results in which viscous effects are treated to 
first and second order. Our main motivation is to understand where hydrodynamics is reliable, how 
its validity and predictive power depends on the initial conditions and which traces are left by dissi- 
pative effects in the data. Fig^shows the time dependence of the temperature for initial conditions 
To = 200 MeV and tq = 0.1 fm/c (left plot) compared to to = 1 fm/c (right plot). Remarkably, 
for early initial times, the first order (Navier-Stokes) equation results in an unphysical increase 
of temperature or total system energy, which we discussed already in the context of Eq. (|34|). In 
contrast, the second order dissipative hydrodynamic equations lead to the monotonous decrease of 
temperature with expansion time, which one expects on general grounds. In general, dissipative 
effects prolong the cooling time. Initializing the system with the same initial temperature To = 200 
MeV at a later time ro = 1 fm/c leads to significantly smaller dissipative effects, see Fig. ^ In 
this case, the Navier-Stokes first order solution provides a quantitatively reasonable estimate for 
the causal second-order Israel-Stewart theory. This raises the question how the applicability of 
hydrodynamics depends on the initialization time ro. 

A. Breakdown of hydrodynamic picture at small To 

The two initial conditions used in Fig^ ar e difficult to compare, since the systems differ signif- 
icantly in entropy density. For the calculations in Fig. we work with initial conditions for which 
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FIG. 2: The freeze-out time Tf (at which the system reaches Tf D = 150 MeV) as a function of initial time 
To for ideal (dashed line) and second order dissipative hydrodynamics (solid line for r/fs = 0.3 and dotted 
line for rj/s = 0.08). The initial temperature T was chosen such that Tq t = const, see text for details. 

Tq tq = const ~ (251 MeV) 2 . For zeroth order ideal hydrodynamics, this ensures that the choice of 
the initialization time tq does not affect the time evolution of the system, see (I32j) . Irrespective of 
To, the system has the same entropy, and thus leads to a final state with the same event multiplicity. 
In this sense, the systems initialized at different tq but Tq to = const are equivalent. 

However, source gradients are larger at earlier times. Therefore, dissipative fluid dynamics will 
reveal stronger deviations from ideal fluid dynamics, if it is initialized at smaller To- This is seen 
in Fig. 12 where we plot the freeze-out time Tf (defined as the time at which the system reaches 
T = 150 MeV) as a function of To- For ideal fluid dynamics, the freeze-out time is unaffected by To- 
In the presence of dissipative effects, however, an earlier initialization provides stronger gradients 
and more time for deviations from ideal hydrodynamics to establish. In general, dissipative fluids 
take longer to cool. The larger the ratio of viscosity over entropy density rj/s, the more pronounced 
is the effect. 

If the gradients in the source are too large, then the gradient expansion underlying dissipative 
hydrodynamics cannot be expected to converge. By estimating uncertainties in 2nd order dissipa- 
tive fluid dynamics, we find that this problem becomes relevant for tq<1 fm/c. At earlier times, 
dissipative hydrodynamics gradually loses its predictive and interpretational power. 

B. Breakdown of hydrodynamic picture at large px 

The validity of a hydrodynamic picture is not only limited to sufficiently late times, but also 
to sufficiently low transverse momenta. To estimate the maximum transverse momentum p™*, 
up to which fluid dynamics may apply, we consider small departures Sf <C 1 of the phase space 
distribution / from its equilibrium value /o- For the gradient expansion underlying dissipative fluid 
dynamics to apply, one requires that local deviations / — /o of the phase space density are small 
compared to the phase space density /. This leads to (see Appendix EJ) 

Sf = ^Tp) U ^ pU « L (41) 
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FIG. 3: The time-dependence of the dissipative strength $/(e + p) (left plot) and the critical transverse 
momentum p™* (right plot). Results are shown for initial conditions To — 300 MeV, To = 1 fm/c, and for 
T]/s = 0.3 (solid line) and rj/s = 0.08 (dashed line), respectively. 



We now consider a particle at point (r, 77, r, (/?), which follows the Bjorken flow, r/ = 
y, and has transverse momentum px- Its four- momentum reads (p l ,p x ,p y ,p z ) = 
(m,T cosh y, pt cos tp p ,PT sin ip p , mr sinh y) , where m T =Vt + m particie depends on the particle mass 
m par ticie- Averaging over the azimuthal angle tp p , we find that (|41j) is equivalent to the condition 

(42) 

Here, $ determines the value of the shear tensor, see (jiffljl . The dissipative strength of the system 
is then quantified by the inverse Reynolds number |48| 

-l $ 

-^Reynolds = £ + ' 

In Figure 01 we plot first a typical time evolution of the inverse Reynolds number ^/(e + p) 
for To = 300 MeV, tq = 1 fm/c and different values rj/s = 0.3 and ry/s = 0.08. Remarkably, the 
value for $/(e + p) shows a maximum at finite time. Therefore, in general, = V4T 2 (e + p)/$ 
assumes its lowest value before freeze-out. Thus, the maximum px at which one can still trust a 
hydrodynamic description may actually be lower then the apparent value obtained at freeze-out. 
This is shown in the second part of FigEl where p™* is plotted as a function of time. As can be 
seen, for moderate values of r]/s ~ 0.3, one obtains a < 0.5 GeV. 

In Fig0J we study how the inverse Reynolds number 3>/(e + p) at the freeze-out temperature 
Tf = 150 MeV depends on the initial conditions tq and To. In general, the dissipative strength 
increases if tq decreases or if Tq decreases. 



C. Hadron spectra at freeze-out 

To calculate hadronic spectra from the phase space distribution / evolved up to the freeze-out 
time Tf , one customarily makes use of the Cooper- Frye freeze-out prescription [43J 

d?N _ fp^ f (44) 



(Pprdy J (2-7T 
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FIG. 4: The dissipative strength $/(e + p) at freeze-out (Tf D = 150 MeV) as a function of initial time To 
(T = 300 MeV fixed, left plot) and as a function of initial temperature T (to = lfm/c fixed, right plot) for 
rj/s — 0.3 (full lines) and ij/s = 0.08 (dotted lines). 



Here, cffi^ denotes the oriented freeze-out volume. For the case of a Bjorken one-dimensional 
expansion, p^dY, 1 * = uit cosh(y — r])rdr]rdrd(j). We now consider a phase space distribution / 
at freeze-out, that deviates locally from equilibrium /q. Adopting the ansatz of Ref. ji^, which 
embeds information about dissipative effects directly in the freeze-out distribution, we write 



/ = fo 1 + 



1 



2T 2 (e+p) 



(45) 



Since the relation (|44jl between the particle spectrum and the phase space distribution is linear, this 
allows us to define the particle spectrum in terms of its ideal part and its dissipative corrections. 
Specializing to the case of a Boltzmann distribution foipuU* 1 ) = 2{N 2 — 1) exp(— p^u^/T), one 



finds for the ideal part [5 



d 2 N° 
d?pTdy 



(2vr)3 /0 



= 2{N 2 C - ^mTTfo^Ki 



V T 



(46) 



The total yield is proportional to the transverse area tt Rq over which we have integrated. For 
numerical calculations, we use Rq = 6 fm. For a Boltzmann distribution and longitudinally boost- 
invariant flow, the dissipative correction to this spectrum reads 



d 2 6N 
d 2 pxdy 



n 



a3 



2T 2 (e 



P^ 



"w 9 



p)J (27T)3^> 

2(W C 2 - l)m T r {o R 2 



4(e + p) 



(2vr)2 

1 (m T \ 
~2 VT) 



Pt\ 2 „ (m T 



V T 



\ T 



(47) 



If r\l s is very small then dissipative corrections will be small. For instance, for the small value 
1)1 1 s = 0.08, one finds a small dissipative strength <&/(e+p) ~ 0.02 at freeze-out, if one uses tq > 0.8 
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P T [GeV] p T [GeV] 

FIG. 5: Particle spectra for pions, kaons and protons from ideal hydrodynamics (open symbols) and viscous 
hydrodynamics (full symbols), for R = 6 fm/c. Left plot: spectra from ideal and viscous hydrodynamics 
($>/(e+p) = 0.02) for the same freeze-out conditions. Right plot: spectra for ideal and viscous hydrodynamics 
(3>/(e + p) = 0.05) with freeze-out of the ideal case adjusted such that slope of viscous calculation is 
reproduced at intermediate pr, (see text for details). 



fm/c and To > 300 MeV for the initial conditions (see FigHJ). This estimate can be reconciled with 
the discussion in Ref. 1431 . if one uses Tf ~ 9 fm/c and Tf D ~ 150 MeV so that Tf Tf ~ 7 instead 
of Tf Tf ~ 1 used in 49]. For this small dissipative strength, we expect only slight modifications 
of hadron spectra. This is confirmed in Fig|SJ 

However, freeze-out time Tf and temperature Tf D are not directly measurable. Rather, they are 
intrinsic parameters of the hydrodynamic model and its possible matching to a hadronic rescatter- 
ing phase. This raises the question to what extent dissipative effects can be accommodated in an 
ideal fluid description of data by varying the freeze-out conditions of the ideal fluid. In Fig. [51 we 
have calculated hadronic spectra iV 1 = N° + 5N for a sizable dissipative strength &/(e + p) = 0.05 
and freeze-out conditions (rf ,Tf ) = (14fm/c, lAOMeV). Remarkably, these spectra can be repro- 
duced to a significant degree by an ideal fluid model, using the same initialization but a different 
decoupling time and temperature (Tf Q , Tf Q ) = (17fm/c, 170MeV). The only clearly visible difference 
between both models in Fig. |S] is a characteristic mass-dependent dip in the low-p± part of the 
spectrum which increases with increasing dissipative strength. This dip is a direct consequence of 
the analytical form (|47|) . and appears to be rather sensitive to small values of the inverse Reynolds 
number. We have observed with curiosity that the fine-binned data available from BRAHMS [5jJ 
for p± > 0.5 GeV for anti-protons maybe hints a similar dip around the lower end of the exper- 
imental reach. We regard this as an illustration that improved particle-identified measurements 
of hadronic transverse momentum spectra in the non-relativistic momentum regime (e.g. smaller 
error bars from PHOBOS [H^]) may provide valuable constraints for dissipative properties of the 
produced dense matter. 



V. CONCLUSIONS 

The predictive and interpretational power of fluid dynamics simulations relies on knowledge 
about the initial conditions, the dissipative properties during evolution, and the modeling of de- 
coupling in the final state, all of which are accompanied by significant uncertainties in phenomeno- 
logical applications. This raises the question to what extent information about fundamental dissi- 
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pative properties of dense QCD matter can be constrained experimentally in heavy ion collisions. 
To address this question, we have compared here results from ideal and dissipative fluid dynamics. 

First, from Boltzmann transport theory (see Section [H] and Appendices), we have established 
that derivations of second order hydrodynamics based on the implementation of d^s^ > 0, (|18|) 
do not provide the complete set of second order terms in a gradient expansion. They discard 
those terms in the gradient expansion which are orthogonal to the shear viscous tensor. We 
then explored in a model study the region of applicability of fluid dynamics. We found that for too 
early initialization times tq<1 fm/c, source gradients are too large for dissipative hydrodynamics to 
converge (see Figs.^and EJ). The inverse Reynolds number $/(e+p) defines a useful p^-dependent 
scale p™* up to which dissipative hydrodynamics can be expected to account for deviations from 
the local equilibrium of a perfect fluid. For our model, we found typically a p^^l GeV. 

Once we had delineated for our model the range of validity of a fluid dynamics description, 
we asked whether data taken within this range can be interpreted unambiguously in terms of the 
strength of dissipative properties of matter. Our results illustrate the challenges of this task. While 
even for a very small dissipative strength, model calculations of transverse hadron spectra show 
a measurable sensitivity to dissipative properties, a very good control over initial conditions and 
freeze-out is required to extract these features from data. In particular, in our model hadronic 
spectra including sizeable dissipative effects could be reproduced rather satisfactorily from perfect 
fluid calculations by modifying the freeze-out temperature and time (see Fig. EJ). 

The fluid dynamics model studied here is based on a strongly simplified transverse geometry 
and does not account for a hadronic scattering phase. Therefore, this setup does not exhibit 
all characteristic features of a realistic heavy ion collision, such as transverse radial and elliptic 
flow. However, the two-step logic used here to assess the range of validity and the interpretational 
power of our model applies also to the more complex, realistic fluid dynamics simulations of heavy 
ion collisions: First, a comparison of ideal and dissipative fluid dynamics results delineates the 
range of validity of a fluid dynamics description. Second, within this limited range, one has to 
establish by model studies to what extent conclusions about the equation of state or dissipative 
properties are independent of the initial conditions and the modeling of the freeze-out. For more 
complex models, there are more measurable quantities (such as elliptic flow), which may help to 
constrain the fluid dynamics. However, there are also more features of the initial conditions, which 
require specification (such as the initial transverse density and flow profile). This increases the 
complexity of the task to disentangle effects of initial conditions from those of the properties of 
matter determining fluid dynamics behavior. We expect that the numerical values for the bounds 
on initialization time and transverse momentum of a reliable fluid description, established here, 
change only weakly by going to more complex fluid simulations. 
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APPENDIX A: DISSIPATIVE HYDRODYNAMICS FROM KINETIC THEORY: 
DERIVATION FOR A MASSLESS BOLTZMANN GAS 

In this Appendix, we derive the equation of motion Q of dissipative hydrodynamics from kinetic 
theory for a massless Boltzmann gas. Our starting point is the Boltzmann transport equation 

p%/(x,t,p) =C(x), (Al) 

with C the collision term. We make use of the notational shorthands ((4 ]) -(|ll |) . We consider the 
first three moments of the Boltzmann equation (|A1|) . using J dio = J 

J cLup^d^f = d^N" = = J dtoC, (A2) 
J doj P »p a d fl f = d^ a = = j duj p a C , (A3) 
J doo p^p^pPd^f = J duj p a pPC . (A4) 

Here, equations ()A2|I and (|A3|) stand for charge conservation and energy-momentum conservation, 
respectively. Equation (|A4|) is the lowest order relation of moments of (|A1|) . which contains dy- 
namic information. For simplicity, we consider in the following an equilibrium distribution /o of 
the Boltzmann type for one degree of freedom, foiu^p 11 ) = exp(— (3u^p^). Also, we employ the 
relaxation time approximation in Eq. (|A4|) . 

C = -p^ f -^ . (A5) 
m 

For all points (x, t, p) in phase space, where we want a hydrodynamic description to apply, we 
assume that departures from local equilibrium are small, 

f = fo(l + Sf). |*/|<1 (A6) 

The correction term 5f can be written in the form 

<5/(x, t, p) = e(x, t) + e A (x, t)p x + e A „(x, t)p x p u . (A7) 

The deviations from the equilibrium energy-momentum tensor Tq U can be likewise written as 

= Tg v + W v + IIA^ + . (A8) 

Here, the deviations from local equilibrium due to shear viscosity, bulk viscosity and heat flow are 
parametrized by H tJ,u , n and respectively. We limit the following discussion to the effects of 
shear viscosity, which is expected to provide for the most significant dissipative effects. It can be 
shown H that in this case we may ignore both the zeroth and first order term in p in Eq. (jA7|) as 
well as concentrate on T <fiu> only. We are left with <5/(x, t,p) = e Al/ (x, t)p x p u , so that the second 
moment of the Boltzmann equation reads 

d^ [ dwp»p <a pP>fo (l + e 75 pV) = / Wp¥/o . (A9) 

To derive from this expression the equation of motion for the non-equilibrium deviation Ii^ v of the 
energy momentum tensor, we establish first the relation between and H^ u . We start from 

< T>» - TP > = J dup<^>(f - f ) 

= e aP {x) J du;p<^>p a p^f 

= e aP (x)I<^>^ , (A10) 
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where 



jMiAtt-im = J dujp^p^ 2 . . .p^f . (All) 

For the Boltzmann distribution /o = ex P ( — P^ ufl /T), the integral IjAllj) can be solved analytically 
for general n |3jj|. In the following, we need the cases n = 4, 5 which read 



I^ a P = a Q vf t u v u a vP + at i^ v u a u P + perm) 

+a 2 (A^A Q/3 + A>" a A ul3 + A^A Q1/ ) , (A12) 

perm. I 

+b 2 (A^A a/? u 1 +perm.\ . (A13) 

The temperature dependent coefficients and hi in (|A12|) and l)A13|) can be calculated e.g. by 
going to the fluid rest-frame, = (1,0). One finds a 2 = 4T 6 /tt 2 and b 2 = 24T 7 /vr 2 . The first of 
these coefficients is related to the pressure p and energy density e as a 2 = T 2 (e + p). Combining 
Eqs. (|A"8j) . (jMOj) and (lA"l2l . we find 

e a/3 (x)/<^ >Q/3 = 2a 2 e a/3 A Q< ^A i/>/3 

= = n<^> . (A14) 

This equation implies u fJi U^ u = and A^IT^ = 0. Therefore, we can set 

e nv = J_ n ^ _ ( A15 ) 

2a 2 



This allows us to rewrite Eq.(jA9J) as 

d ( p<a[)> + jM <a/3>7^ = _ fl^£ J <a/3> 7 5 e (Alg) 
V / r n 

The three terms entering this equation can be rewritten with the help of the identities derived 
above, 

e lS uJ <aP> ^ St = 26 2 e 7(5 A 7<Q A /3>5 = ^n <Q/3> , (A17) 

d M j dup^p <a pP> /o = -I^ <a ^d^ = < V Q ?/ > , (A18) 

I^ <aP> ^d^ 5 = 2b 2 (AT <Q A /3>5 ^d M e 7(5 + 2-uW <a A /3> %e 7(5 ) . (A19) 



Further simplifications arise from u^n 75 = -Ii 1& d^v? and A^rF" 5 = e^A^IF 5 = 0. We 
then decompose Ii^ u d^u y into a totally symmetric component oc B^iy , and a totally antisymmetric 
expression proportional to the vorticity tensor uj u ^ 

IP%« 7 = ir^uv + ft^), (A20) 
= A^ a A^d [p u a] , (A21) 
^ = A^A^d (0 u a) . (A22) 



At this point, we follow [36| in assuming the 'rigidity of flow'. This means that one assumes 
the equilibrium flow u M to be shear free, and thus terms proportional to 0^ v to vanish. As a 
consequence, the term e 7( j J dujp^p <a p /3> p 1 p s d fJi fo can be neglected. Then, Ea. (|A16|) reads 

— < V"/ > +— f A^A^DU^g + 2ir 5 ( Q a/c )N ) + b 2 U af3 Da^ = - — —U a ^. (A23) 
T a 2 V / m 02 
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For the third term on the left hand side, Da^ 1 oc De oc d^u^ oc 9fi, which vanishes under the 
assumption of rigid flow. We note that in Ref. pvR IsH ]. this term is kept. The contractions 
u v d^T llv = and A va dnT ,a ' = for an energy momentum tensor of the form T^ v = Tq V + Ii^ v 
lead to the equations of motion Q and (J2J), respectively. To complete the derivation of ©, we 
require that in the limit rn — > 0, the first order dissipative Navier-Stokes equations are recovered 
from (jA2.3|) . This leads to the identification 

4^ ^ * (A24) 
The equation of motion (j3J) then follows from ()A23|) . 



APPENDIX B: DISSIPATIVE HYDRODYNAMICS: DERIVATION FOR A MASSLESS 
BOSE-EINSTEIN GAS AND MOMENTUM-DEPENDENT RELAXATION TIME 

In this Appendix we derive the second order hydrodynamic equation (J3J for the case that the 
relaxation time rn is momentum-dependent, rn = p^u^f, f = const. This ansatz is of interest 
since it often enters calculations of the shear viscosity (see Refs. H E3, H and sec. 8.3 of the 
textbook |Hj|). In contrast to Appendix we start here from the equilibrium distribution for a 
massless Bose-Einstein gas /q = [exp(/3u M p^) — 1] , having in mind QCD with gluonic degrees of 
freedom. This implies that the departures from equilibrium are expressed by [c.f. Ea. l|A6|) ] 

<5/(x, t, p) = (1 + / )eA,(x, t)p X p u . (Bl) 

The relation between e^ u and is expressed as in Eq. (|A15|) . except that the coefficient a 2 
is replaced by 02 = 4C(5)T 6 /7r 2 , with ((n) the Riemann ("-function. This implies that 0,2 = 
^|yT 2 (e + p) when expressed in terms of energy density and pressure of the free Bose-Einstein gas. 
This modification, compared to the Boltzmann gas used in App[^J is due to the replacement of 
the integrals Jw«2-/% ^ 

/W2-Mn = J duJ p^l p H2 _ . . p M"/ (l + / Q ) . (B2) 

The second moment of the Boltzmann equation (|Alj) reads now 

^ (><a/3> + eiS P<*P>l8\ = _hS_f<aP> 7 S ? ( g 3) 



where the right hand side depends on an order four tensor only, in contrast to the order five tensor 
entering Ea. ()A16|) . Consequently, for a Bose-Einstein gas we find 

d,I^ = -^<W>, (B4) 
e 75 I <Q/3>7<5 = n 75 A 7<a A^ >5 = 2n a/3 , (B5) 
p<ap> 7 S d ^ s = 2 " h (A^A^^d^s + 2uW <Q A /3><5 d M e 7(5 ) , (B6) 



where 



24 

6 2 = -C(6)T 7 . (B7) 

7T Z 



Together with the assumption of the "rigidity of flow", i.e. setting Q^ v = 0, Eq.(jB3j) becomes 

Zh. < V <V > +^ (A^A^DU^s + 2U S ^) = -~^IT*f (B8) 
T a,2 v 1 ra 2 
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Thus, we recover Eq. (J3J) upon identifying 



V <-> 7F T ' r n <~> — T, (B9 
T a 2 



which leads to the ratio 



b 2 T _ 3 C(4)C(6) ^ in24 3 

In conclusion, the dissipative hydrodynamic equations for a massless Boltzmann and Bose-Einstein 
gas differ, but only on the order of a few percent. 



APPENDIX C: PERTURB ATIVE CALCULATION OF THE RELAXATION TIME 

To complete the discussion on viscosity and relaxation time, we calculate in QCD the relaxation 
time f = Tuip^u^) -1 to leading order in the strong coupling constant g. We follow the line 
of argument of Refs.0, [j^ . l55| for a gluonic system with 2(iV 2 — 1) = 16 degrees of freedom. 
The collision term C, entering this calculation, is then dominated by 2 <-> 2 gluon scattering 
processes. The corresponding effective matrix element M gg is given in the hard-thermal loop 
(HTL) approximation, including Landau damping To simplify the calculation, a small time- 
independent velocity u x (y) is considered which only varies with the space coordinate y 0, 

p M #/(x, t, p) ~ -/ (1 + /o)^/3^^ • (CI) 

Following analogous steps as in 0, EH)j t can be written in terms of the ratio 



1 2(iV2 _ i) 



f 16 



f[ J du Jl (27r) 4 5\p 1 +p 2 -p 3 -p 4 )\M gg \2(pf p y i ) 2 fo(p 1 

-i=l 



x/o(P2) (1 + /o(Ps)) (1 + /o(P4))] / (/ Mp x P y ?fo(l + /o)) • (C2) 



Here, |M 99 | 2 denotes the scattering matrix element squared, summed (averaged) over spins and 
color degrees of freedom in the final (initial) state. The integrals are evaluated e.g. in [55J], leading 
to 

1 (TV 2 - 1W 5 T 2 (/ , . ; 2x 

?- 120 c(5) (W 1 -' 4 ^ '- (C3) 

From Ea. ()B9|) . the estimate for rj is derived 

when taking the color degrees of freedom into account. We note that 0,2 is proportional to a factor 
2 (TV 2 — 1). This expression agrees also with the value quoted in 
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